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ABSTRACT 

We investigate Gamma-ray Burst (GRB) internal shocks with moderate magne- 
tization, with the magnetization parameter a ranging from 0.001 to 10. Possible mag- 
netic dissipation in the stripped magnetized shells is also taken into account through 
introducing a parameter k {0 < k ^ 1), which is the ratio of the electric field strength 
of the downstream and the upstream. By solving the general MHD jump conditions, 
we show that the dynamic evolution of the shock with magnetic dissipation is different 
from the familiar one obtained in the ideal MHD limit. As long as the relative velocity 
between the two magnetized shells is larger than the corresponding Alfven velocities 
in both shells, strong internal shocks still exist for a ^ 1, which can effectively tap 
kinetic energy into radiation. However, in the ideal MHD limit (fc — 1), the upstream 
magnetic energy can not be converted into the downstream thermal energy so that 
the GRB radiation efficiency is low. This is likely inconsistent with the current GRB 
data. With magnetic dissipation, e.g., k ^ 0.5 the range of k is constrained given a 
particular upstream-downstream Lorentz factor 721 and a magnetization parameter 
a), a significant fraction of the upstream magnetic energy can be converted into the 
prompt 7— ray emission. At the typical internal shock radius, the characteristic syn- 
chrotron emission frequency in the magnetic dissipation dominated case is however 
too large (oc cr^) compared with the data if cr ^ 1. On the other hand, as long as the 
ordered magnetic field is stronger than or at least comparable with the random one 
generated in the internal shocks, a net linear polarization ^ 30% results. In view of 
the possible high degree of linear polarization of GRB 021206 and the identification of 
a possible highly magnetized flow in GRB 990123 and GRB 021211, we suggest that 
a mildly magnetized internal shock model (0.01 < cr < 1) with moderate magnetic 
dissipation is a good candidate to explain the GRB prompt emission data. 

Key words: Gamma-rays: bursts-Magnetic fields-Magnetohydrodynamics (MHD); 
shock waves; relativity 



1 INTRODUCTION 

Tremendous advances to understand the Gamma-ray bursts 
(GRBs), one of the greatest enigmas in high energy as- 
trophysics have been achieved in the past seven years (see 
Meszaros 2002; Cheng & Lu 2001; Zhang & Meszaros 2004 
for reviews). However, the nature of the GRB central engine 
is still unclear. In the conventional fireball model, a GRB 
is powered by the collisions of non-magnetized shells with 
variable Lorentz factors, i.e., the internal shocks (Paczynski 
& Xu 1994; Rees & Meszaros 1994; Kobayashi, Piran & 
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Sari 1997; Pilla & Loeb 1998; Daigne & Mochkovitch 1998; 
Guetta, Spada & Waxman 2001). The magnetic field in- 
volved in the synchrotron radiation model for the prompt 
7— ray emission is usually interpreted as being generated 
in the internal shocks, and is randomly oriented in the 
shock plane with small coherence scale, so that there is 
no net polarization expected in the prompt 7— ray emission 
(Medeveder & Loeb 1999). Recently two pieces of indepen- 
dent evidence suggest that the GRB central engine might 
be strongly magnetized. First, the detection of the very high 
linear polarization of GRB 021206 (Coburn & Boggs 2003, 
but see Rutledge & Fox 2004) suggests that the magnetic 
field involved in the synchrotron radiation could be glob- 
ally ordered (e.g. Coburn & Boggs 2003; Lyutikov, Pariev 
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& Blandford 2003; Granot 2003; Granot & Konigl 2003), al- 
though some alternative explanations such as the Compton 
scattering model (e.g., Shaviv & Dar 1995; Lazzati et al. 
2000; Eichler & Levinson 2003) and the narrow jet model 
(Waxman 2003) remain. Second, analysis of the two well- 
studied optical flashes from GRB 990123 and GRB 021211 
reveal that the magnetic fields in the reverse shock region 
are stronger than that in the forward shock region, so that 
the GRB outflows are likely magnetized (Fan et al. 2002; 
Zhang. Kobayashi & Moszaros 2003; Kumar & Panaitcscu 
2003) . The current models involving ordered magnetic fields 
for GRBs invoke a Poynting flux dominated outflow (e.g. 
Usov 1994; Thompson 1994; Mcszaros & Rccs 1997; Spruit, 
Daigne & Drcnkliahn 2001), in which the ratio of the elec- 
tromagnetic energy flux to the kinetic energy flux of the 
baryons (i.e., a) is of order 100 or more, and the GRB 
prompt emission is envisaged to be due to some less familiar 
magnetic dissipation process. In reality, a GRB event likely 
involves a "hot component" as invoked in the traditional fire- 
ball model (e.g. due to neutrino annihilations), whose inter- 
play with the "cold" Poynting flux component would allow 
the cr value to vary in a wide range (e.g., Rees & Meszaros 
1994; Zhang & Mcszaros 2002). On the other hand, numer- 
ical simulations and statistic analysis of GRBs suggest that 
the internal shock model is preferred (e.g., Kobayashi et al. 
1997; Lloyd-Ronning, Petrosian & Mallozzi 2000; Guetta et 
al. 2001; Amati et al. 2002; Zhang & Meszaros 2002; Wei & 
Gao 2003). Motivated by these facts, in this paper, we inves- 
tigate the GRB internal shocks with magnetization. Spruit 
et al. (2001) have discussed this topic briefly by taking the 
energy and momentum conservation of two magnetized shell 
collision. More detailed treatments axe needed. 

This paper is structured as follows. We first discuss the 
MHD 90° shock jump conditions with magnetic dissipation, 
both analytically (§2.1) and numerically (§2.2). We then (§3) 
discuss the fast ejecta - slow cjocta interaction, in particu- 
lar for an ejecta with moderate magnetization (i.e., a < 10), 
and calculate the prompt synchrotron emission from such 
moderately magnetized internal shocks. Our results are sum- 
marized in §4. 



2 INVESTIGATION OF THE RELATIVISTIC 
MHD 90° SHOCKS WITH MAGNETIC 
DISSIPATION 

It is generally believed that a GRB involves a rapidly ro- 
tating central engine. If the magnetic fields from the engine 
are frozen in the expanding shell, the radial magnetic field 
component decreases more rapidly with radius (oc r~^) than 
the toroidal field component (oc r~^). At the internal shock 
radius, the frozen-in field is likely dominated by the toroidal 
component, so that the field lines are essentially perpendic- 
ular to the shock normal direction, i.e. one has a 90° shock. 

2.1 The general Jump Conditions 

A rigorous analytical treatment of the MHD 90" shock jump 
conditions has been presented by Zhang & Kobayashi (2004, 
hereafter ZK04) recently. Similar to ZK04, here we present 
a rigorous analytical solution for the 90° shock jump condi- 
tion with magnetic dissipation. The main difference between 



this work and ZK04 is the following. In their ideal MHD 
(i.e., without magnetic dissipation) limit, the electric field 
in the shock frame is continuous across the shock. This is no 
longer true in the presence of magnetic energy dissipation, 
which may result, e.g., from magnetic rcconncction in the 
shock front (e.g. Levinson & van Putten 1997; Lyubarsky 
2003). In this work, following the treatment of Lyubarsky 
(2003) on the termination shock in a stripped pulsar wind, 
we introduce a parameter 

k = f32,B2s/l3i,Bu, O^k^l (1) 

to describe the potentially important but poorly under- 
stood magnetic dissipation process, where pis and /32s (-Bis 
and B2s) are the velocities (magnetic field strength) of the 
upstream and downstream regions measured in the shock 
frame, respectively. In principle, k is constrained by the 
stripped part of the Poynting flux, which is assumed as a 
free parameter ranging from to 1 in the current work. 

Now, following ZK04, we consider a relativistic shock 
that propagates into a magnetized ejecta. In the following 
analysis, the unshocked region (upstream) is denoted as the 
region 1, the shocked region (downstream) is denoted as the 
region 2, and the shock itself is denoted as "s"'^. Hereafter 
Qij denotes the value of the quantity Q in the region "i" in 
the rest frame of "j" , and Qi denotes the value of the quan- 
tity Q in the region "i" in its own rest frame. For example, 
721 is the Lorentz factor of regions 2 relative to region 1, 
(3is is the velocity (in unit of the speed of light c) of re- 
gion 1 relative to the shock front, B2s {E2s) is the magnetic 
(electric) field strength of the region 2 measured in the rest 
frame of the shock, while Bi is the comoving magnetic field 
strength in the region 1, etc. In the presence of the mag- 
netic energy dissipation, the familiar relativistic 90° shock 
Rankine-Hugoniot relations (Hoffmann, De & Teller 1950; 
Kennel & Coroniti 1984, hereafter KC84) take the general 
form 

mUls = n2'ti2s , (2) 

-Els = PibBib; E2b = /32s-B2s , (3) 
7isMi + = 72sM2 + , (4) 

UlUls + + = W2^t2s + -22- + ELIEL (5) 

where P denotes the dimensionless velocity, 7 = (1-/3^)"^''^ 
denotes the Lorentz factor, and u = fS'f denotes the radial 
four velocity. Hereafter, n, e, p — (F— l)e denote the number 
density, internal energy and thermal pressure, respectively, 
and F is the adiabatic index. The enthalpy is nrripC^ + e+p, 
and the specific enthalpy can be written as 

^ = ^^c^ + (6) 

where rup is the proton mass and c is the speed of light. It 
is convenient to define a parameter 

"^i = ~A 9" 5 (i) 

VRfii^i 47rni/ii7;g 

Notice that such a notation system is only vahd for §2 and the 
Appendix. When discussing the GRB problem, i.e. the fast shell- 
slow shell interaction (§3), we introduce different meanings for 
the subscript numbers. 
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to denote the degree of magnetization in each region. The 
magnetization parameter in the upstream region (cri) is a 
more fundamental parameter, since it characterizes the mag- 
netization of the flow itself. We therefore define 

a = a^= ^ ^ . (8) 

In our problem, we are interested in a "cold" upstream flow, 
i.e., ei = pi = 0, so that /ii = rripC? . This is the only 
assumption made in the derivation. With equation Q, the 
thermal Lorentz factor of the downstream particles is 

+ (9) 

712 mpC^ r 72s P2s 

where 725(721, f^, fc) is a function of 721, o, and k, and can 
be solved once 721, o and k are known. The equation gov- 
erning 725(721, fT", fc) reads (see Appendix for derivation, see 
also Lyubarsky 2003) 



1 fr-1) 1 1 

7i5Ui5(l + -)[/3is - /32s - ^ , ^ 1 + ^ 



= k 



(2-r)n?, (f-l)Mi, 



2r uo 



Fcr U2s 



(10) 



where ^is = (/92s+/32i)/(l+/32i/32s), 71s = 72s72i(l+/32i/32s) 
and Mis = /3is7is (e-g-, ZK04). 

Now, the compressive ratio can be derived directly from 
equation ||5J, i.e. 



^^^^^ [nL(72i,a) + l]^/^ ^^2^_^^i,2^(^^) 



112s (721, (j) 



712 _ «ls (721,0-) 

ni M2s (721,0-) 

In the downstream region, the total pressure includes 
the contribution from the thermal pressure P2 = (f — l)e2 
and the magnetic pressure pb,2 = B2s/87r7|s. The ratio be- 
tween the magnetic pressure to the thermal one is: 

\P2sJ ^87r(f-l)e27|s 
k^a Mis 



P2 



62 



(12) 



2(r - 1) M2s 'n2?npc2 

where equations ^ and Q have been used. 

For simplicity, in §2, we take T = 4/3. For 721 S> 1, 
7is ~ 72s 721(1 -I- /32s) > 721 > 1- Equation can be 
solved analytically (see also Lyubarsky 2003) 



where y = -^-2- . We then have 



n2 721(1 + /32s) 7 + X + + 14x + 

/^2s 1 + X + \/l + 14x + ' 



(13) 



(14) 



62 



^(o + l)(l+/32s)(l--f ) 
4 P2s 

721(1 + a) 
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(7 + X + V1 + 14X + X^ 

6x 1 



i + x+v/i + i4x + : 



(15) 



2X 



Pb,2 ^ 

P2 ~ (/32s - X) 



12X 



[1 - 5x + + 14X + X' ] 



(16) 



Since e2/n2m-p(? > 0, equation 1151 hints that /32s > X- Foi' 
o- 3> 1, 721 S> 1 and fc = 1, X ~ 1 — ^, equations p3l - ll6^ 
are reduced to 



P2 



i 712 
2(7 ' Til 



272 



62 



7l27npC^ 



3 Pb,2 ^ . .,_s 
7721, -4(7,(17) 

4 P2 



so that 712/(4721 + 3)711 1/2, 62/(721 — l)7i27npC^ a; 3/4. 
All these are well consistent with the numerical results pre- 
sented in ZK04 (see also our Figure0(b) and (c)). This hints 
that strong shocks still exist in the high a regime, as argued 
by ZK04. On the other hand, without magnetic dissipation 
(k = 1), the upstream magnetic energy can not be converted 
to thermal energy and radiation. Since Pb,2/P2 « 4o- ^ 1, 
in the high-cr regime, only a tiny amount of total energy 
(kinetic plus magnetic) can be used for electron synchrotron 
radiation. The GRB radiation efficiency in the A; = 1, a ^ 1 
model is therefore too low to interpret the data, which is 
typically above 40% (Lloyd- Ronning & Zhang 2004). 

For 721 1 and x ^ 1, equations 113I - I16II are reduced 

to 



P2 



l + 4x 



7l2 
Til 



4721 , 



62 



n277lpC-^ 



i 721(0 + l)(l-2x). 



Pb,2 
V2 



6x. 



(18) 



For (J <C 1, equations II8II are reduced to /32s 
2 

i.pC 



1/3, 7I2/7I1 

4721, 62/7i27npC^ « 721 and Pb,2/p2 « 60 ~ 0. All these are 
consistent with the familiar results presented in Blandford 
& McKee (1976). As shown in equations 11711 . for fc = 1, 
721 > 1 and (7 » 1, Pb,2/P2 ~ 4(7. So, roughly speaking, 
for fc = 1 and 721 ^ 1, Pb,2/p2 is linearly proportional to a, 
which is the case (see the thin dotted line shown in Figure 
Hd)). 

When writing down eas. l2l5ll . one has already assumed 
that a pair of shocks form as the collision happens. In order 
to have a shock, the relative velocity of shells (the corre- 
sponding Lorentz factor is 741 following the standard con- 
vention for shock jump condition analysis) should be faster 
than the Alfven velocity {v\ — c^J the corresponding 

Lorentz factor 7a = \/l + o). So the condition is 741 > 7a. 
If this condition is not satisfied, there will be no shock wave 
and the two shells would simply bounce off of each other 
elastically. 

If a shock form, in the rest frame of the upstream the 
velocity of the shock front should be faster than ma, i.e., 
7is = 72s 721(1 -I- /32s/32i) > 7A. In the ideal MHD limit, 
this is always the case. However, in the presence of mag- 
netic dissipation, 71s > 7a is not satisfied for an arbitrary 
fc (0 ^ fc ^ 1) value. This in turn imposes a strict con- 
straint on the possible choices of fc. This can be understood 
as follows. In the presence of magnetic dissipation, espe- 
cially for fc ^ 1, /32s ~ 1/3 or even smaller, 71s > 7a yields 
721 > ^W"^ J bearing in mind that 721 is always smaller 
than 741. Therefore for a given small value of fc, 721 has to 
be larger than a certain value in order to have a physical 
solution of the shock. In other words, a certain value of 721 
and o, fc is constrained within a certain range. When fc = 1 
(no dissipation), the shock can happen in the whole 721 — 0- 
plane as long as 741 > Vl + o is satisfied. 
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Figure 1. The variations of four parameters, i.e., U2b, 
e2/n2mpC^, n2/ni and Pb 2/P2, as functions of 721 (sub-figure 
a, b, c and d, respectively). The thin hnes are for fc = 1, i.e., the 
ideal MHD case. The thick lines are for the case with significant 
magnetic dissipation, i.e. k = 0.5. The solid, dashed, dotted, dash- 
dotted and dash-dot-dotted lines are for cr = 0.01, 0.1, 1, 10, 100, 
respectively. Similar to ZK04, the parameter e2 /n2mpC^ (thermal 
Lorentz factor in the shocked, downstream region) is normalized 
to (721 — 1), and the parameter n2/ni (compressive ratio) is nor- 
malized to (4721 -1-3), both being the values expected in the a = 
case. For the cr = 1, 10, 100 cases, the e2/n2mpC^ values are too 
large to fit into the scale, and we have multiplied the values by 
a factor 1/3, 1/20, and 1/200, respectively (see sub-figure c). For 
a given fe < 1 and cr, the condition 71s > 7a is not satisfied un- 
less 712 is larger than a critical value. So some lines do not cover 
the whole horizontal axis scale. Similar situations also happen in 
Figure 2. 



2.2 Numerical Investigations 

For general cases with arbitrary 721, a and fc, equation 
iiave to be solved numerically, and equations IIH and 
II2II can be calculated correspondingly. 

Figure Q and |21 siiow tiie variations of four parame- 
ters, i.e., U2b, e2/n2mpC^, 712/711 and Pb,2/P2 as functions 
of 721 (Figure and a (Figure respectively. The thick 
lines are for fc = 0.5 (the case with significant magnetic 
dissipation), and the thin lines are for the k = 1 (the 
ideal MHD limit case). In Figure the thick solid, dashed, 
dotted, dash-dotted and dash-dot-dotted lines are for a = 
0.01,0.1,1,10,100, respectively. The thin solid and dotted 
lines are for a — 0.01, 100, respectively. The parameter 
e2ln2rap(? (thermal Lorentz factor in the shocked, down- 
stream region) is normalized to (721 — 1), and the parameter 
n2/n\ (compressive ratio) is normalized to (4721 -I- 3), both 
being the values expected in the cr = case (Blandford & 
McKee 1976). For the a — 1, 10, 100 cases, the e2/n2mpC?' 
values are too large to fit into the scale, and we have mul- 



Figure 2. The variations of four parameters, i.e., U2b, 
e2/n2mpC^, 712/711 and Pb 2/P2, as functions of cr (sub-figure a, 
b, c and d, respectively). The thin solid and the thin dotted line 
represent fc = 1 and 721 = 1.5, 100, respectively. The thick solid, 
dashed and dotted lines are for k = 0.5 and 721 = 1.5,3.0, 100, 
respectively. Similar to Figure^ the parameter €2 / n2mpC^ is nor- 
malized to (721 — 1), and the parameter 712/711 is normalized to 
(4721+3). 



tiplied the values by a factor 1/3, 1/20, and 1/200, respec- 
tively. In Figure |5| the thick solid, dashed and dotted lines 
are for fc = 0.5, 721 — 1.5, 3, 100, respectively, and the 
thin solid and dotted lines are for fc = 1, 721 = 1.5, 100, 
respectively. Similar to Figure the parameter e2/n2mpC? 
is normalized to (721 — 1), and the parameter 712/711 is nor- 
malized to (4721 -f- 3). 

In Figure when 721 > 10, nearly all the (normalized) 
parameters are insensitive to 721 (see also ZK04 for the case 
of fc = 1), but are sensitive to fc and/or a (see also Fig- 
ure Inland Figure |HJl. For instance, as shown in Figure Q^c), 
in the absence of magnetic dissipation, 62/(721 — V)n2mp(? 
is of order unity for a G [0.01,100]. But for fc = 0.5, 
62/(721 — l)Ti2mpC^ <x {a + 1) increases linearly with cr (see 
Figure |5Jc)). This is because at higher a values, more and 
more magnetic energy is converted to thermal heat while 
the total number of leptons is decreasing with a, so that 
the energy per electron increases rapidly with o (see also 
Zhang & Meszaros 2002). For fc = 1, ct > 1 and 721 > 1, 
It 2s — \f^- But for fc = 0.5, the increase of U2s with o is only 
slightly (see Figure0^a),|5Ja) for detail). A similar situation 
is evident in FigureQ^d) and|2Id). For fc = 1, Pb,2/P2 ~ 4o- 
increases linearly with o. However, for fc = 0.5, the resulting 
Pb,2/p2 is always of order unity for a ^ 1, which hints that 
significant part of upstream magnetic energy can be con- 
verted into the downstream thermal energy, no matter how 
large the cr is. 
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Figure 3. The variations of two parameters, i.e., pj, 2/P2 and 
e2/n2mpC^, as functions of k. The solid line and dotted line are for 
721 = 2.5, 100, respectively. The thin lines represent a = 10, the 
thick lines represent (7 = 1. The parameter e2/n2mpc'^ (thermal 
Lorentz factor in the shocked, downstream region) is normalized 
to (721 — 1), the values expected in the cr = case. 



Figure 4. The condition 71s > 7a poses a constraint on the 
range of fc as a function of 721 and cr. The solid, dotted, dashed 
and dash-dotted lines are for cr = 0.3, 1, 10, 100 respectively. 
Above each line is the physically allowed region. 



To see the impact of magnetic dissipation on the dy- 
namic evolution of the shock more clearly, in Figure El we 
plot the variables e2/n2m-p(P' and Ph,2/p2 as functions of k. 
The solid and the dotted lines are for 721 = 2.5, 100, re- 
spectively. The thin lines represent cr = 10 (For a = 10, 
(t/(1 + (7) « 1, so the thin solid line and dotted line nearly 
overlap), the thick lines represent cr = 1. The parameter 
e2/n2'7ipC^ has been normalized to (721 — 1), the values ex- 
pected in the cr = case. The results presented in Figure 
are consistent with our equations 1181 . e.g., for k < 0.5 and 
cr > 1, e2/n2mpC^ ^ {1 — 2x)(cr + 1)721 » (721 — 1) and 
Pb,2/P2 ~ 6x oc fc^. 

As mentioned before, the condition 71s > ja imposes a 
strict constraint on k for a given value of 721 and cr. Here 
we discuss it in more detail. The condition 71s = 725721 (1 + 
/32s/32i) > 7a is plotted in Figure 0] where the solid, dotted, 
dashed and dash-dotted lines are for cr = 0.3, 1, 10, 100 
respectively. Above each line defined by a particular cr and 
an arbitrary k, a dissipative shock is allowed. Given a par- 
ticular cr value (i.e. for a particular line), one can constrain 
the fc value range given a 721 value. For a small enough 721, 
each 721 value corresponds to a minimum k value, so that 
the shock can not be too dissipative. As a decreases, the 
restriction is progressively weaker. 



3 INTERNAL SHOCKS POWERED BY THE 
COLLISION BETWEEN TWO 
MAGNETIZED SHELLS 

3.1 The internal shocks 

We now use the results obtained in §2 to study internal 
shocks. Consider a faster, trailing, ultra-relativistic, magne- 
tized shell (Ff) hits a slower, leading, magnetized shell (Fs), 
where Ff and Fs represent the Lorentz factor of the fast and 



slow shells (measured by the observer), respectively. The 
corresponding velocities for the two shells are /3f and /3s, 
respectively. Upon collision, two shocks form, i.e. a reverse 
shock (RS) propagating into the fast shell and a forward 
shock (FS) expanding into the slow shell. The shocks in- 
crease the thermal and magnetic densities of both shells. 
There are four regions in this system, i.e. the unshocked 
slow shell (1), the shocked slow shell (2), the shocked fast 
shell (3), the unshocked fast shell (4). A contact discontinu- 
ity separates the shocked fast shell material and the shocked 
slow shell material. In the following analysis, velocities Pr^ 
and their corresponding Lorentz factors Fi = (1 — P'ri)~^^^ 
(Fi = Fs and F4 = Ff) are measured in the observer frame. 
Thermodynamic quantities, e.g., rii, pi, ei (particle number 
density, thermal pressure, thermal energy density) are mea- 
sured in the rest frame of the fluid, so are the magnetic 
pressure and magnetic energy density, i.e. pe.i and eB,i. In 
this work, we assume that these two shells are cold, i.e., the 
specific enthalpy /ii — fi4 = nripC? (see equation (|SJ for the 
definition) . 

The equation that governs the FS takes the form (with 
equation llOt 'l 



7itti(l + 

cri 

,^2 (2 - f2) uj 
= K2 o 

2F2 ui. 



/32- 

(f2 



(r2 



1) 1 



F2 

1) m 



U272 



+ 



F2Cri U2 ■ 



(19) 



where cri = i3i/[47r7f nimpC^], B\ (-B2) denote the FS frame 
magnetic field strength of region 1 (2); r2 is the adiabatic 
index of the region 2; = f^^' 
Similarly, for the RS we have 



74^^4(1 H )[/34 

0-4 



(r3 



U373 2 



6 



2r3 



(ra - 1) M4 

f 30-4 ^3 ' 



(20) 



where 0-4 = i3|/[47r7|n4mpC ], ,63 (-B4) denote the RS frame 
magnetic field strength of region 3 (4); T3 is the adiabatic 
index of the region 3, ^3 = Here 71 and 72 (73 and 

74) are the forward (reverse) shock frame Lorentz factor 
of the fluids in the region 1 and 2 (3 and 4), respectively, 
and u? = 7j^ — 1 (j = 1 — 4) are the 4-velocities of the 
fluids, /3j = Mj/7j^ The 7j's can be parameterized by the 
Fj's as follows, i.e., for FfajFrs 3> 1 and Fj ^ 1, one has 
271 ~ (Ffs/F, + F,/Ffe), 272 ~ (Fts/Fa + F2/Ffe), 273 « 
(r,e/r3+r3/F„) and 274 ~ (F.e/Ff + rf/r,,). Here Ffa and 
Frs are the Lorentz factors of the FS and RS measured in the 
observer frame, respectively. These equations in turn suggest 
that F3 « (74 — M4)(73 + W3)Ff and r2 « (72 — M2)(7i+iii)Fs. 
Therefore the equality of the velocities along the contact 
discontinuity (r2 = F3) yields 



Ff ~ (74 + U4)i'y3 — U3)(72 — '"2)(7l + Ui)rs 



(21) 



The equality of the total pressure (the sum of the thermal 
pressure and the magnetic pressure) along the contact dis- 
continuity (P2,tot = -P3,tot) yields 



{%^S(i + '-4(i-fcif|)) 



+ 



fe^}n3M4.(22) 



Equations 119I22II are the basic formulae for the following 
calculations. 

In this work, it has been assumed that two shocks form 
in a collision between two magnetized "shells". The condi- 
tion for this to happen is that the relative velocity of the 
two shells exceeds the va- If the relative Lorentz factor be- 
tween the two shells is 741, the condition can be written as 
741 > VI + o"- 
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Figure 5. The bulk and thermal Lorentz factors and the mag- 
netic equipartition parameter as functions of a. Ff/Fs = 20 is 
adopted throughout. The thick lines represent k = 0.5, the thin 
lines represent A; = 1, i.e., the ideal MHD limit, (a) The solid 
lines represent the bulk Lorentz factor of the shocked region, i.e. 
r2 (= r3), normalized to the value for cr = 0. For k = 0.5 and 
1, the resulting F2's are similar, so that the thin/thick solid lines 
overlap. The dotted/dashed lines represent the thermal Lorentz 
factors of the shocked baryons in FS/RS regions, respectively, 
again normalized to the values for cr = 0; (b) £3 1 the downstream 
magnetic energy density normalized to the sum of the downstream 
thermal energy density and the magnetic one. The solid/dotted 
lines correspond to the eb obtained in the FS/RS regions, respec- 
tively. 



3.2 Numerical Results 

The problem is complicated and numerical calculations are 
needed. For simplicity, we assume that the magnetization 
parameters a, the dissipation parameter k, the rest masses, 
and the widths of the shells (measured by the observer) are 
the same for both shells, so that cti = 0-4 = a, k2 = ka = k, 
Pi = P4 = and / = Ui/ni = Fs/Ff < 1 (which im- 
plies that the RS is stronger than the FS). As a numerical 
example, we assume a ranges from 10~^ to 10, and take 
Fs = 50, Ff/Fs = 1// = 20, i.e., 741 ~ 10. We extend the 
a range up to 10 in the following calculation, so that the 
shock form condition 741 > \/l -|- a is always satisfied. Our 
calculations therefore satisfy the shock form condition. The 
dissipation parameter k is taken for two indicative values, 
i.e. 1 and 0.5. Another important parameter involved is the 
adiabatic index of the shocked material. Here we treat it 
self-consistently: we define the thermal Lorentz factor of the 
downstream baryons as 7th, i' = e-^i jn-^im^c? {i'=2, 3). For 
7th, i' > 1 we take Fj/ = 4/3 since both electrons and protons 
are relativistic. Otherwise, protons are only sub-relativistic 
although electrons are relativistic. In such a case, we have 

^ Notice that the notations here are different from those in §2 in 
that we have dropped out the subscript "s" here for simplicity. 



f;/ = 13/9. The equations 1191221 are then solved numer- 
ically. With the resulting 7j (j = 1 — 4), we can calculate 
r2, 7th, i' as a function of a and k directly. The numerical 
results are shown in Figure (K^a). Since we are mainly inter- 
ested in the novel features introduced by the magnetic fields, 
we normalize the relevant parameters with their correspond- 
ing value for cr = 0. For k — 1 and a ~ 10~^ — 10, the results 
change only slightly with respect to those for cr = 0, which 
implies that strong internal shocks still exist in the high cr 
ZK04 argued recently (see the thin lines plotted in 
Figure 1^ a) for detail). For k = 0.5, the results are quite 
different. For example, 7th, i' oc a when k = 0.5, which can 
be very high if a ^ 1 (see the thick lines plotted in Fig- 
ure [KJ a) for detail). The main reason is that for a constant 
k, the total dissipated energy essentially remains the same 
for different cr values while the number of leptons decreases 
sharply as oc cr~^, so that the energy per lepton increases as 
(X a. 

The downstream magnetic field is amplified effectively. 
For convenience, we define the ratio of magnetic energy den- 
sity to the sum of the thermal energy density and the mag- 
netic one as 

eB,i' = ; 

eB,i' + Si' 
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Figure 6. The net linear polarization degree, the radiation effi- 
ciency, and the typical synchrotron frequency as functions of a. 
The thick lines represent k = 0.5, and the thin lines represent 
k = I, i.e., the ideal MHD Umit. Tf/Fs = 20 is adopted through- 
out, (a) The solid line represents the net linear polarization of the 
emission from region 3, as a function of cr (where £3,0 = 0.01 is 
adopted); The dotted line represents the efficiency of converting 
the upstream total energy into the prompt synchrotron emission 
frequency as a function of cr, as normalized to the value for cr = 0. 
(b) The observed typical synchrotron radiation of the shocked 
electrons in region 2 (the dotted line) and 3 (the solid line). Other 
parameters adopted here are Sc = 0.5, R = 10^'^ cm, £30 = 0.01, 
p = 2.5, and Lq = lO^^ergs. The cosmological redshift correction 
effect is not taken into account. 



^ ^"i' (23) 

^ + Ti7[^(l + -(l-fc^f^))-l] 

Please note that unless a <^ 1, it is much different from the 
familiar notation eb.o, the fraction of the shock-generated 
magnetic energy density to the total shock-dissipated energy 
density in the non-magnetized fireball model. Here i" = 1, 4, 
which corresponds to i' — 2, 3, respectively. The numerical 
results are shown in Figure |SJb). For k = 1, the magnetic 
energy is not dissipated, so the downstream magnetic energy 
is strong. For instance, for cr ^ 3 x 10"'^, the downstream 
eB,i' ^ 10~^ , which is strong enough to match what is needed 
in the internal shock synchrotron model of GRBs (Guetta et 
al. 2001). For cr ^ 1, Eb,!' reaches a asymptotic value ~ 1.0, 
since the downstream magnetic energy density dominates 
the total ones. However, for k — 0.5, i.e., significant part of 
magnetic energy has been dissipated, the resulting magnetic 
energy is much weaker than that of fc = 1. For cr ^ 1, the 
resulting Eb,!' reaches a asymptotic value ~ 0.4. 

If our line of sight to the ejecta is not very near the edge 
of the cone, due to the beaming effect, the viewed outflow is 
essentially axis-symmetric. As a result, the net polarization 
contributed by the random magnetic field is nearly zero. 
On the other hand, the existence of the ordered magnetic 



field component in the ejecta likely results in a net linear 
polarization (e.g. Granot 2003) . The observed net linear po- 
larization can be expressed as (Granot & Konigl 2003; Fan 
et al. 2004) 

Hnet- 0.60^3^, (24) 

where is the ratio of the ordered magnetic energy den- 
sity to the random one. In the current case, for the re- 
gion 3 (Relatively speaking, the radiation comes from the 
region 2 is weaker and softer, which contributes little to 
the 7— ray emission. So we mainly focus on the region 3), 
— £3,3/(1 — £8,3)53,0. We take the typical value of the 
random field equipartition parameter as £b,o ~ lO"'^ (e.g. 
Guetta et al. 2001). The net linear polarization degree is 
calculated in Figure [HJ a) (the solid lines). We can see that 
only cr > 3.0 X 10~^ is required to produce Ilnot > 30% for 
the k = 1 case, while for k — 0.5 one requires a > 0.02. This 
more demanding requirement is simply due to that much 
of the downstream magnetic energy has been dissipated for 
k = 0.5. 

Assuming that the shocked electrons in the region i' 
have a power law distribution in energy, i.e. dn/d'y^y oc 

(7c,i' > 7c,i',m) withp > 2, wehave7e,i',m « £e7th,i'[(P- 
2)/(p — l)](mp/mo). As usual, we assume that the inter- 
nal shocks take place at a radius R ~ 10^'^cm and the 
wind luminosity is Lq = lO^'^ergs s~^. We can then cal- 
culate the typical synchrotron radiation frequency v^ny ~ 
£^i'7o i' ^^2/'2TvmcC. Both the ordered magnetic field com- 
ponent and the random magnetic component are taken into 
account. 

Another important parameter involved is the conver- 
sion efficiency of the internal shock, i.e., the fraction of 
the total upstream energy (kinetic plus magnetic) that is 
converted into the downstream thermal energy. This quan- 
tity is directly related to the GRB radiation efficiency 
through the parameter te, i.e. the energy which is trans- 
ported to electrons can be radiated effectively in the fast- 
cooling regime (which is justified in the GRB prompt emis- 
sion phase). Since the emission of the region 3 is much 
stronger and harder than that of the region 2 (see|Sfb)), 
here we mainly consider the former component. The down- 
stream thermal energy per baryon (in the observer frame) 
is eth,3 ~ 7th,3r2mpC^, and the total number of baryons in- 
volved is Nt = LoSt/r{{l + a)mpC^, where i?tot = LoSt is 
the total upstream energy (in the observer frame), and the 
factor 1/ ( 1 + cr) represents the fraction of the kinetic energy 
to the total energy. The conversion efficiency can be then 
written as e = eth,3A^6/-Btotee7th,3r2/rf (1 -I- cr). This effi- 
ciency is plotted against cr in Figure |nj a) (dotted lines) for 
both fc = 1 (thin dotted) and k — 0.5 (thick dotted). Since 
we are mainly interested in the novel features introduced 
by the magnetic fields, the efficiency has been normalized 
to the corresponding value for ct = 0. We can see that e 
generally decreases with a. For k = 1, the decrease is dras- 
tic so that the radiative efficiency for a ^ 1 is too low to 
interpret the GRB data (Panaitescu & Kumar 2001; Lloyd- 
Ronning & Zhang 2004). For k — 0.5, one can still retain a 
high efficiency in the high-cr regime. 

In Figure|SJb) we plot the typical synchrotron frequency 
as a function of cr for both the RS and the FS (each calcu- 
lated for both fc = 1 and k = 0.5). We find that for the typi- 
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cal parameters adopted here, the forward shock emission has 
both a lower frequency and a lower luminosity than the RS. 
For fc = 1, the RS emission frequency is about 2 x lO^^Hz, 
which is insensitive to a and matches the current obser- 
vation. However, as discussed above, the a ^ 1 regime is 
disfavored due to the low efficiency involved. For k = 0.5, 
the downstream magnetic energy can be converted into the 
prompt 7— ray emission effectively, but for a ^ 1, the result- 
ing RS emission frequency i/^i oc is much harder than the 
BATSE band (Fig|^b); see also Zhang & Meszaros 2002). 



4 DISCUSSIONS 

In the standard fireball model of GRBs, the prompt 7— ray 
emission is believed to be powered by internal shocks (e.g., 
Meszaros 2002). That model is successful to interpret the 
GRB variability (e.g., Kobayashi et al. 1997) and some em- 
pirical relations (e.g., Zhang & Meszaros 2002). However, 
if the magnetic field involved in the internal shock region 
is random (e.g. generated within the shocked region due to 
plasma instability) , it is difficult to account for the observed 
high linear polarization of GRB021206 (Waxman 2003; Gra- 
not 2003). If the involved shells are magnetized, even if the 
magnetization parameter is only mild (e.g. a > 10~^ with- 
out magnetic dissipation and a > 0.02 for substantial mag- 
netic dissipation), the amplified downstream ordered mag- 
netic field is strong enough to dominate the random com- 
ponent and generate a significant degree of linear polariza- 
tion (see Figure|SJa)). For a lower magnetization parameter, 
e.g., a <g 10~^, no observable net high linear polarization is 
produced unless some other geometry effects are taken into 
account. 

Through introducing a "magnetic dissipation parame- 
ter", k (Lyubarsky 2003), the MHD jump condition with 
magnetization is solved, and the parameterized expres- 
sions and numerical calculations for the compressive ratio 
(nd/riu), the magnetic pressure - to - thermal pressure ra- 
tio (j3b,d/Pd), and the downstream mean random Lorentz 
factor (ed/nampC^) (where the subscript u, d represent the 
upstream and downstream respectively) are presented for 
various 721, cr and k values (See §2). In this work, we treat 
a and k as independent parameters and investigate the de- 
pendence of the solution on both parameters. As expected, 
we found that the results obtained for shocks with magnetic 
energy dissipation is very different from the familiar one 
obtained in the ideal MHD limit. The introduction of the 
k parameter manifests our ignorance of the poorly known 
magnetic dissipation process, and its value is poorly con- 
strained. In reality, k may be correlated with a and 721 . For 
example, in the a ^ limit one has fc ^ 1. Also in Fig. 
|l|we have shown that the lower limit of fc is jointly deter- 
mined by 721 and a. However, lacking a theory for magnetic 
dissipation, we simply treat fc as a free parameter as long 
as it satisfies the condition shown in Fig. |1| With the gen- 
eral jump conditions, the GRB internal shocks with moder- 
ate magnetization (10~^ < a < 10) have been calculated. 
Various considerations/constraints allow us to narrow down 
the a range for a possible GRB model. We show that for 
fc = 1, i.e., the ideal MHD limit, strong internal shocks still 
exist in the high a case (a 2> 1). However, in the ct 2> 1 
regime, the upstream magnetic energy can not be converted 



into the downstream internal energy effectively, resulting in 
a very low radiation efficiency inconsistent with the current 
GRBs data (Panaitescu & Kumar 2001; Lloyd-Ronning & 
Zhang 2004). For a significant magnetic dissipation case (e.g. 
fc ^ 0.5), the upstream magnetic energy can be converted 
into the prompt 7— ray emission effectively for an arbitrary 
cr, but for (7^1, the observed frequency (oc a^) is much 
harder than what we observe giving a typical internal shock 
radius. We therefore disfavor a a value much greater than 
unity. We note that too high a typical radiation frequency 
is a common feature for any high-cr model, and the prob- 
lem may be remedied by considering a possible pair cascade 
in the magnetic dissipation region (e.g. Zhang & Meszaros 
2002; the pair emission in the internal shocks can be found 
in Fan & Wei 2004 and Li & Song 2004). Developing a de- 
tailed pair-dominated model in the high-a regime is however 
beyond the scope of the present work. 

At the low-cr side, if the claimed high linear polarization 
in GRB 021206 is true (Coburn & Boggs 2003; cf. Rutledge 
& Fox 2004), within the synchrotron model, the required 
magnetization parameter is cr > 10~^ for fc = 1 and a > 0.02 
for fc = 0.5 in order to give rise to a ^ 30% linear polariza- 
tion. Modelling early afterglows for GRB 990123 and GRB 
021211 generally requires a magnetized flow (Fan et al. 2002; 
Zhang et al. 2003; Kumar & Panaitescu 2003), with cr > 0.1 
(ZK04) for these two bursts. Considering more general cases, 
we favor a mildly magnetized fireball model (10~^ < cr ^ 1), 
first suggested in Rees & Meszaros (1994). A traditional 
problem of the internal shock model is its low radiation ef- 
ficiency (e.g. Panaitescu, Spada & Meszaros 1999; Kumar 
1999). In the presence of a toroidal magnetic field, the effi- 
ciency of the internal shock is even lower in the ideal MHD 
limit. In view that the magnetic dissipation process (which 
is naturally expected for a moderate cr value) can help to 
solve this problem (see FigurelSJa), the thick dashed line for 
detail), we suggest that a mildly magnetized internal shock 
model with moderate magnetic dissipation is a good candi- 
date to explain the current GRB prompt emission data. 

The mild magnetization (e.g. cr ~ (10~^ — 1)) preferred 
in this paper is likely a natural outcome of a realistic cen- 
tral engine. In a GRB event, a rapidly rotating magnetar- 
type (either black hole - torus system or neutron star) cen- 
tral engine with a surface magnetic field ~ 10^^ G likely 
launches a Poynting flux flow with an isotropic luminos- 
ity ~ 10^" — lO^^ergs s^^. The cataclysmic event also in- 
volves a hot flreball component due to processes such as 
neutrino annihilation, with a typical isotropic luminosity of 
~ 10^^ — lO^'^ ergs s~^. The latter energy component may be 
orders of magnitude stronger than or is at least comparable 
with the former, so that the picture recommended here is 
justifled. 
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APPENDIX A: DERIVATION OF EQUATION 
(10) 

Equations l|U and @ can be rearranged as follows (with 
niuis = n2M2s, (Al) 

7lsUlsWl = 72s«2siy2, (A2) 

wiul, + {pi + B'(/8tt) = i«2ML + (p2 + B|/87r), (A3) 

where wi = n^i + Bf /4tt = nimp(? + jrTYPi + Bf /Atv, Hi 

(i=l,2) are measured in the comoving frame, Fi are the adi- 
abatic index of region i. With the definition of Wi, we have 

fi - 1, 



Pi — — — — -{wi — mmpf? — Bf /in). 
Equation 1A2II then reads 

r 2 ^ ri-1 , fi-1 2 ^ 2-fi Bl 

Wi(WisH ) nimpC H — 

Fi Fi 2ri 47r 



(A4) 



2 , 2 — f 2 



^2,r2 — 1. r2 — 1 2," ^ z /■ \ r\ 

^W2(u2s^ ) n2TnpC H — . (A5) 

r2 r2 2r2 47r 

On the other hand, equation <A3ll yields W2 = 
7isitisWi/72s'it2s. Substituting this into equation 1A5II . we 
have 



7isUisi(;i(/3is — /32s + ^ 



Fi-l 1 Ta-l 1 



Pi 71bU1b To 72a"2B 



-nimpC -j^n2mpC + ^jr-— 



For pi <^ ninipC^, which we are interested in here, one has 
wi ^ nimpC^(l + a). Considering n2/ni = U\a/u2a and 
[B2/Bi_f = fc2[■ule/^t2s]^ we have 



7isUi,(l + cr)(/3is-/32s + 4 



Fi-l 1 Tj-l 1 



Ti 71s"ls r2 723"23 

_ r2-i tiiB _|_ 2^^2.1,2^ t£iB^2^ _ 2-f 1 „ 



After some simple algebra, we finally have 

7lsUls(l + o){f3xs ~ /32s 



r2-l 1 'I I a_ 
fa 72a "2b ' 2 



r2-l "la I 2-r2 7.2/Uia \2 



(A7) 



(A8) 



which is equation llOII in §2. 



